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CHAPTER I 


INTRODUCTION 

There are two alternative approaches to the study of the laws of crack 
propagation and toughness characterization. The fracture mechanics 
approach is concerned with the stress intensity factor K or energy release 
rate J which are derived using the concepts of continuum mechanics. The 
material science approach concentrates on knowledge of the hierarchy of 
defects, their development, and interactions; that is, this approach empha- 
sizes the microraechanisms of fracture processes. In the current work one 
of the possible ways of unifying these two approaches based on the CL 
theory is discussed. 

The analysis of the critical energy release rate G c shows that this 

parameter is history dependent. This is also supported by our findings. 

In addition, observations of the kinetics of crack growth show nonmonotonic 

crack advance under monotonic changes of energy release rate J (or stress 

intensity factor K) . This implies the existence of intrinsic properties of 

the material which are not reflected by the G , K , J and K parameters. In 

c c 

many studies on the morphology of regions around the crack [1-6], the 
complex processes of the generation and development of microdefects on 
different hierarchial levels are observed. Perhaps parameters like K and 
yield stress are not sufficient to describe these fracture processes. 
For example, we observe entirely different shapes of the damage zones for 
different steels under similar loading conditions. 

The recently proposed crack layer (CL) theory [7,8,9] considers the 
crack together with the surrounding cloud of defects as one system which 
has several degrees of freedom. In the CL an active zone, where the 
nucleation and development of defects occur, can be distinguished. The 
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active zone (denoted in the literature as a process zone) and the wake zone 
are both distinct In the CL. 

The motion of the active zone Is decomposed into translation, rotation, 
and deformation. The generalized forces associated with the above mentioned 
degrees of freedom are derived within the framework of the thermodynamics 
of irreversible processes. These generalized forces are represented by 
linear functions of path-independent integrals , L, M [10] and by 

integral characteristics of damage R^, R^-jj > R 0 * Consequently the CL 

theory defines the relationship between the parameters of fracture 
mechanics and the characteristics of microstructural changes which are the 
subject of material science. Experimental study of the CL theory 
predictions and formulation of a suitable damage parameter are the ultimate 
goals of this work. 

The toughness characteristic is represented as the product of the 

specific enthalpy of damage and integral cross section damage R^. The 

former (enthalpy) is a candidate for being the material constant according 

to the CL theory. The latter (R^) is history-dependent and is responsible 

for the widely observed changes of G . 

c 

Introduction of the R^ parameter and its evolution are based on the 
analysis of the morphology of defects. The calculation of J or K is based 
on stress analysis and macroscopical measurements, for this reason comple- 
mentary macro- and micro-experiments were carried out simultaneously. 
These experimental procedures are described in Chapter II. The results of 
the observations are presented in Chapter III. Chapter IV contains the 
analysis of the results, comparisons with the theoretical predictions, and 
interpretation of the data. 

Dr. A. R. Rosenfield of Battel le-Columbus Laboratories made many helpful 
suggestions . 


2 



CHAPTER II 


DESCRIPTION OF EXPERIMENTAL PROCEDURES 

This study establishes the relationship of material toughness and the 
macroscopic process of crack growth with microscopic changes. The research 
described herein consists of both macro- and microscopic studies and the 
micro-macro relationship. 

A. The macroscopic studies included inquiries into the characteriza- 
tion of material toughness and the kinetics of crack propagation. 

1 . Toughness Characterization Test 

Samples with cracks of specified length, grown under sinusoidal 
tension-tension loading conditions, were subjected to the standard test for 
fracture toughness evaluation. This procedure uses three samples to obtain 
one data point. The AISI 302, 305 stainless , and 1070 carbon steels were 
used in this investigation. Data on fracture toughness parameters usually 
show large scatter [11]. In order to obtain statistically representative 
data, an ensemble of identical samples was prepared. The geometry of the 
samples is 150mm in length and 20mm in width. The thickness for each part- 
icular material was fixed. For different materials the thickness varied 
from 0.15mm to 0.5mm. A sixty degree notch was cut in the middle of a long 
side of each specimen. The ensemble included eight groups of 10 samples 
each. In each group of samples, cracks were grown to a specified length 

using sinusoidal tension-tension testing conditions (R «* a . /a = 1/3, 

min max 

-2 

average stress O = 150 MNm , frequency v» 50 Hz, room temperature T = 
2 ®°C). All the samples were then subjected to the standard test for frac- 
ture toughness evaluation. The maximum pulling force was recorded automa- 
tically. The critical crack length was measured on the fracture surface 
where the boundary between fatigue crack growth and unstable failure could 
be clearly identified. 
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2. Test for Kinetics of Crack. Layer Propagation 

The study of CL propagation was conducted on AISI 302 and 304 stainless 
steels and 1070 carbon steel. The geometry of the samples is the same as 
described above. A 20 KN capacity electrohydraulic closed loop MTS fatigue 
machine was used to conduct the constant amplitude, sinusoidal waveform, 
tension-tension cyclic loading. The samples were electrolytically polished 
prior to fatigue testing in Jacquet reagent [12] to reduce surface stresses 
and to obtain a mirrored surface for microscopic studies. The crack 
advance was registered on the attached scale using a travelling optical 
microscope. This microscope was attached to the MTS machine (see Figure 1) 
to control and register crack and surrounding damage evolution. 



FIGURE 1 

MTS fatigue machine with attached travelling 
optical microscope. 
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Microscopic studies consisted of direct observation of the CL 


propagation; analysis of dislocation density, and study of discontinuity 
surfaces. The damage zone surrounding the crack was observed through 
metallographlc and SEM microscopes at magnifications of 50X, 100X, 250X, 
500X, 1000X, 10,000X and 20,000X in order to identify quantitatively the 
primary elements of damage. A general picture of the CL was constructed 
from fragments of both side and top views. 

A profile of dislocation density was obtained by a systematic mapping 
of the hardness measurements taken in the region of the crack tip [5, 13, 
14] . Figure 2 shows the grid of the actual measurement points in the 
vicinity of the crack tip. 



Picture of the actual grid consisting of microhardness 
dents on the surface of the fatigue crack specimen in the 
region of the crack tip. The spacing between the centers of the 
indents in x and y directions is 100 

Discontinuities around the crack, which can be identified at 250X and 
higher magnifications were measured as follows. A grid consisting of two 


5 





orthogonal sets of parallel lines was superimposed on the picture taken at 
500X magnification (Figure 3). Discontinuities had a prefered orienta- 
tion. Therefore, the method of biased sampling was employed to obtain 
statistics of discontinuities [15]. The density of discontinuities is cal- 
culated on the basis of the number of intersections of the testing line 
with the traces of discontinuity surfaces [2,16,17]. 



FIGURE 3 

Grid for statistical analysis of discontinuous 
surfaces . 



CHAPTER III 


EXPERIMENTAL RESULTS 


A. Phenomenological Study of Toughness Characterization 

The goal of the current study is to examine the loading history-depen- 
dence of G . We consider the simplest loading history: the constant 

c 

amplitude sinusoidal tension-tension loading. For these loading conditions 

the maximum local stresses in the vicinity of the growing crack tip 

monotonically increase with the crack length. Therefore, the J Integral 

alone is sufficient to characterize the particular loading history. We use 

the value of J at the end of the fatigue stage 'J ' as the parameter to 

reflect the history of loading. Figure 4 shows the dependence of the 

critical energy release rate G on J, for three different materials. This 

c h 

dependence was obtained in accordance with experimental procedures outlined 
in the previous section. The points represented by symbols in Figure 4 are 
the results of three independent tests. 

The critical energy release rate was calculated using the following 
expression: 


G - o 2 7T £ F 2 U /b)E _1 (1) 

c c c c 


where o is the failure stress, F(£/b) is a geometrical factor calculated 
c 

as , 


2b ^ t\ 1 
l 8 2b 


0.752 + 2.02 (l/b) 


+ 0.37 (1-sin ~) 


cos 


ttS, 

2b 


FU/b) 


( 2 ) 





Critical Energy Release Rate , G P (kJm ’ 2 




and b Is the width of the specimen [18]. 

The value of J. Is calculated according to the following formula, 
h 

J, = a 2 tt? F° ( i’ /b)E _1 
h max 

where a x is the maximum applied stress during the fatigue cycle, i is the 
length of the interrupted fatigue crack. 

The statistical analysis of fracture toughness was conducted for AISI 
304 cold rolled, full hard austenitic stainless steel. Figure 5 shows the 
results of 80 experiments on fracture toughness for that material. 
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B. Kinetics of CL Propagation 


The rate of crack growth for AISI 304 cold rolled, full hard steel 

under the previously described conditions was plotted as a function of J in 

Figure 6. After the initial acceleration the crack stopped growing at 

2 

approximately 8 kJ/m , the value of J (see Figure 6). Such a decelaration 
is an often observed phenomenon. When the decelaration of the crack growth 
occurs, the damage zone begins to expand. In the literature this 
phenomenon is called "crack arrest”. The damage expansion (formation of an 
£-zone) associated with the "crack arrest" is shown in Figure 7. 

The general views of the well developed crack layers in 301, 304 and 
1070 AISI steels are shown in Figures 8, 9 and 10. These pictures ex- 
plicitly show that crack propagation is accompanied by a surrounding layer 
of intense 'damage'. 

Macroscopic configurations of the CL appear to vary among various 

materials tested under identical loading conditions. The only parameter 

2 

with length dimension in the conventional fracture mechanics is (K/a ) . 

y 

-2 -2 

The yield stresses a for the tested materials are 780 MNm , 920 MNm and 

y 

-2 

590 MNm . The experimentally measured non-dimensional parameter r /{, 

P 

takes the following values at the 6mm crack length: 0.0341, 0.0476, 0.095. 

2 

At the same time the theoretical characteristic size (K/o ) /£ for the 


y 

considered conditions equal correspondingly 1.6, 1.18 and 1.22. Therefore, 
K and a are not sufficient to describe the observed phenomena. This 

y 

discrepancy calls for additional parameter that would reflect differences 
in the CL growth process on the microscopic level. Such a parameter, 
damage density P, will be introduced in the next Chapter. 
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Energy Release Rate JikJm -2 ) 


FIGURE 6 

Crack growth rate vs. energy release rate in 304 AISI 
stainless steel. The decelaration part of the curve 
corresponds to the formation of £-zone. 
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FIGURE 8 


Fatigue crack in 0.2 mm thick specimen of 301 AISI stainless steel. 
Frequency of loading is 50 Hz, max stress is 25% of the yield stress 
(R = 0.2). Note the damage layer around and ahead of the main crack 
Crack length - 9mm. 



FIGURE 9 


Fatigue Crack in 0.15mm thick specimen of 301 stainless steel. 
Frequency of loading is 50 Hz, max stress is 25% of the yield 
stress (R = 0.14). Crack length - 7mm. 



CRACK TIP 


i - 1 

oi 



Fatigue crack in 0.5mm thick specimen of 1070 AISI steel. 
Frequency of the loading is 50 Hz, maximum stress is 25% 
of the yield stress (R = 0.2). Note the damage layer 
around and ahead of the main crack. Crack length - 7mm. 



C . Morphology of Crack Layer 

If the damage is viewed on a progressively finer scale a hierarchy of 
defects can be visualized. At low magnification nonhomogeneous deformation 
of continuous media is observed; no discrete defects could be identified 
(Figure 11). In Figure 12 the traces of discontinuous surfaces, which 
correspond to slip band extrusions and intrusions, are clearly observed in 
304 AISI stainless steel. Further magnifications of the object shown in 
Figure 12, which is seen as a single line, show that it represents on 
average 20 discontinuous surfaces (Figures 13a and 13b). If we are to 
continue to view the damage at the progressively higher resolutions the 
following elements could be observed: clusters of dislocations, single 
disloctions, subgrain boundary precipitates, point and lattice defects. In 
this paper two types of defects were considered and the energy dissipation 
associated with these were estimated. 

The first type of defect to be investigated is the distribution of dis- 
locations around the crack tip since dislocations are commonly associated 
with the plastic deformation. In Figure (14) the map of dislocation density 

, around the fatigue crack based on the measurements of microhardness is 
disi 

shown. In order to plot this map the following relationship between the 
microhardness (DPH) and dislocation density P ^ ^ , obtained on the basis of 
[19,20,21] was used: 



where and C 2 are experimentally obtained constants. A three-dimensional 
picture of dislocation density around fatigue crack is shown in Figure 15. 
The discontinuity surfaces represent the second type of defects under 
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consideration. A statistical analysis of the discontinuities was done ac- 
cording to the methodology described in Section II B. 


The discontinuity surface density P disl *- s introduced by the formula 

[16,17] 


P 


disl 



(4) 


where P is the number of point intersections per unit length of test line. 

1 5 2 3 

The density of discontinuity surfaces changes from 10 to 10 mm /mm at the 
observable regions. A map of the discontinuity distribution based on this 
methodology is shown in Figure 16. 


18 




FIGURE 1! 


Morphology of crack layer at different magnifications: 
(a) General view of the fatigue crack at low magnifi- 
cation; (b) The crack tip region enclosed in the box on 
11a. Extensive damage is seen around and in front of 
the crack tip. The features of an element of "damage" 
become clearer under higher resolution (see Figure 12). 






FIGURE 12 


Picture of the crack tip region enclosed in the box shown in Figure lib. 
Extensive damage consisting of slip band intrusions and extrusions is 
clearly observed. The arrows indicate position of the wake of the crack. 

The crack itself is out of focus because the focusing was aimed at the 
slip bands. Pictures at this magnification were used for statistical eval- 
uation of discontinuity density p around the crack. The finer structure 
of the slip bands enclosed in the box become clearer under higher resolution 
(See Figure 13a). 













FIGURE 13 


SEM pictures of the crack tip region taken at 10,000 and 20,000 
magnifications.. The crack tip is shown by the arrow in Figure a. 

In Figure b the region enclosed in the box in Figure a is magnified 
The traces of the discontinuity surfaces are clearly visible. 


FIGURE 14 

Contours of equal levels of dislocation density around the 

crack tip in 304 AISI stainless steel. Contour A corresponds 

7 -2 

to the dislocation density 10 cm . B*A+d, C**A+2d, D*A+3d , 

7 -2 

Q»A+19d, where d>*5.10 cm . Symbol "+" indicates the crack 



FIGURE 15 

cm 

Three dimensional picture of dislocation density p[— 3] 
around the fatigue crack in 304 AISI stainless steel. 
Sample No. 174. Crack length - 7mm. 
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FIGURES 16 

.mm 2 , 

The countours of equal levels of discontinuity density pi — ;J 
around the fatigue crack m 304 AISI stainless steel. 

Symbol "+" indicates the crack tip. 
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Region of the quasistatic crack growth in AISI 304 
stainless steel. The arrows indicate the typical 
microcrack propagation in the direction of the main 
crack. 


Enlarged region of the area enclosed in the box in 
Figure a. The surface has a typical brittle frac- 
ture appearance. 



CHAPTER IV 


DISCUSSION 


Our results show that the critical energy release rate G c depends on 
the methods employed in the preparation of the specimen for the fracture 
toughness test. This agrees with the results of other workers. A 
statistical analysis of the data from Figure 5 is presented on Figure 19. 



Energy Release Rate 

FIGURE 19 

Critical energy release rate as a function of the energy release 
rate J^. Solid line represents the average values of G c# Dashed lines 
indicate the 63% confidence zone. This plot is result the of statisti- 
cal evaluation of the data shown in Figure 4. The value of 2.5kJm""^ 
corresponds to minimal crack length. 
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A mathematical expectation is shown by the solid line and the 63% confi- 
dence zone is shown by the dotted lines. 

As illustrated by this plot, the average energy release rate G in- 

c 

creases with the increase of J, . It is important to note that the variance 

n 

of also increases with the increase of J^. The latter is an important 
factor for engineering design applications. 

In this study we observed nonmonotonic crack advance under monotonic 
changes of energy release rate J. There is an apparent alternation between 
the crack advance and damage growth processes. This implies the existence 
of intrinsic properties of the fracture process which are not reflected by 
the J integral. 

The monotonic crack growth as well as the dependence of G and its 

c 

variance on the history of the process call for an explanation in terms of 
the morphologic changes accompanying the crack propagation. One possible 
explanation is given in the CL theory. 

A. Concept of Crack Layer 

CL is represented schematically in Figure 20. 



Sketch of crack layer. 


28 



The zone surrounding the crack where P > P Q (where P Q is the reference 

level of damage) is called crack layer. The active and wake zones are 

distinct in the following. The active zone where the damage growth takes 
• • 

place (P>P Q , p>o where P is a rate of damage growth) is bound by the lead- 
ing and trailing edges ^ , T ^ ^ . xhe wa k e ZO ne can be described as the 
trace of the active zone movement and characterized by p>p Q » P ■ o . We 
are considering only two elementary movements of the active zone: transla- 

tion as a rigid body and isotropic expansion (the spread of damage), 


v fc (x) - \( x °) + e(x-x^) 


( 5 ) 


o 

where v, is a rate of the translation of the center of the ative zone x , 
k 

x-x° is a radius vector for arbitrary point of the active zone, e is a 
rate of isotropic expansion. (We look upon the distribution of defects as a 
mass distribution and the center of the active zone is found as a center of 
mass . ) 

The measure of isotropic expansion 'e' is given by the logarithmic 
measure 


e » In d/d 


( 6 ) 


where d is a characteristic size of the active zone and d is that size in 

o 

the initial configuration. 

For rectilinear CL growth 

v^ » £ - rate of crack length growth 

v 2 - 0. 

Then according to the CL theory the rates of extension (l) and expansion (e) 
are given as 


01 , 


J(£)d£ 


l = 


yR 1 - J^Ji) 


( 7 ) 
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and 


Si 


£ 

J(£)d£ 


e = 


_o 

yR - M 
o 


( 8 ) 


where J, and M are the well known path independent integrals: 


J - | | » 1 


J, = 


'V "i - a ij n j u i,i ldr; 


(9) 


(*) 


M = 6 


k£ 


w 


l f £ Vl ' °ij :) J u i.k Jt d dr 1 


the resistance moment 

*i\ 


u) 


p *n^dT , 


( 10 ) 


( 11 ) 


stands for a measure of resistance of the crack extension, and 


R 


o 


J pdA, 

A 


( 12 ) 


is a measure of resistance for crack expansion. 

Above f, n, , o and u stand for strain energy density, k-th component 

0 K. lj X 

of the unit normal vector with respect to a contour T, stress tensor compo- 
nents, and displacement vector component correspondingly. stands for 

the derivative with respect to x^ , and 6 is a Kroneker's delta. From the 
Equation (7) the rate of crack growth becomes uncontrolled (l»® ) when 


approaches R 


1* 


Thus, conventional G could be identified with R : 

c 1 

M £*) = G . (13) 

c 


Therefore, history dependency of G^ corresponds to the evolution of the 
resistance moment R^. The density damage can be represented as 
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(14) 


P “ <5(x 2 > + p(x 2 > 

* 

where stands for core of damage associated with the crack edges and 

their immediate vicinity (boundary layer), the Dirac's delta func- 

* 

tion which indicates that the core of damage R is concentrated on the 
crack trajectory (x 2 = 0) and P(x^,x 2 ) stands for the measured damage 
distribution outside of core of damage. 

Integrating both sides of Equation (14) over the trailing edge T(t) 

we 

find: 


where R^ is the variable part of the resistance moment. 

In the regions where fracture surface has the same appearance one would 

* 

expect the core of damage to be the same. 

The f ractographic analysis of the fracture surfaces shows no changes in 

the fracture appearance in the region of the quasi-static crack growth. 

* 

Therefore, we assumed R^ to be a constant. The value of this constant is 

associated with microcracks indicated on Figure 18. Study of the sizes, 

orientation and density of these microcracks for quantitative evaluation of 
* 


R^ is under way. 

According to the equations (13) and (15) the critical value of G can 

c 

, , . . . . _ associated with the core of damage R* and 

be resolved into two parts: Gc s 1 


G c associated with the damage dissemination around the crack, 

G = G* + G (16) 

c c c 

The damage density G c introduced above represents either dislocations 
or discontinuity surfaces, or any other type of defect or their superposi- 
tion. 
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The value of G is measured in the conventional fracture toughness 
c 

evaluation (macroscopic) test. On the other hand, according to the CL 

theory, G can be evaluated in terms of the integral damage R obtained 
c l 

from microscopic analysis the and specific enthalpy > of damage under 
consideration. The comparison of the two values of G^ will reveal the 
weight of each particular defect in the fracture process. Following this 
argument two types of defects, namely dislocations and discontinuity 
surfaces, are considered. 

B. Evaluation of Energy Stored in the Dislocation Network 

Energy E., , of one dislocation could be estimated as 10-20 nJ/m [22]. 
disl 

When the dislocation network is generated the energy associated with 
dislocations per unit crack increment is 


E = E R 
disl 1 


Figure (14) shows the map of the dislocation distribution in the vicinity 
of the crack tip. 

Resistance moment R^ for the map shown in Figure (14) is R^ ~ 10^ m 
Then the energy stored in the dislocation network 
E~10~ 8 J/m * 10 10 m" 1 = 10 2 J/m 2 

constitutes 0.1% from the experimentally measured critical energy release 

2 

rate G c = 120 kJ/m . It is worth noting that the more appropriate proce- 
dure would involve consideration of the entire dislocation network since it 

has been established that there is a definite stage of dislocation density 

12 -2 

saturation for a given steel (10 cm ). But, even this will not change 
much the estimated value of E. 

12 

If we consider the case when the maximum dislocation density (10 

_2 

cm ) is everywhere in the CL, the full energy of dislocation network will 
not be higher than 1% of the G^ value. Therefore, the generation of dis- 
locations can explain neither the value of G nor the evolution of G with 
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crack advance. We conclude that such a simplistic dislocation approach 
cannot account for actual values of critical energy release rate G c « 

C. Evaluation of Energy Associated With the Discontinuity Surfaces 

The integral cross-section damage R^ for 304 AISI steel was experimen- 
tally obtained according to the methodology and formula (Eq. 11) given 
above. Using the method of least squares the functional relationship 

between R and J, is found and plotted in Figure 21 as a solid line, stars 
1 n 

represent experimental data. The experimental data on fracture toughness 

(Figure 5) are replotted in Figure 22 as G vs. R. . 

c l 

The specific energy of a discontinuity surface is presented by the slope 
of the solid line (which is in its turn obtained by the least square 
method) . 

Assuming the average thickness of the slip lines (layers) to be i ■ 5oX 

i 3 ~ 

and mass density d - 8 g/cm , one can obtain the total mass Rj of damaged 
material per square centimeter of crack surface: 

>\j -2 . -2 

R^ = R^6d = 10 g'cm 

't n, 9 -1 

The ratio G c /R = Y represents the specific energy of damage (4.10 .erg.g ) 

which is the same order of magnitude as the heat of fus ion of steel 
9 -1 

(2.9.10 erg.g ) [23] . This result indicates the relationship between fracture 
and a phase transition. The relationship between the processes of fracture and 
melting (and/or dissociation) was noted first by Born [24] and then was 
experimentally studied by Born and Furth [25, 26] . 

Following Born’s idea the physical interpretatin of fracture and 
plastic deformation as an anisotropic localized phase transitions is 
discussed in [27]. 
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FIGURE 21 

Resistance moment R vs. energy release rate. Stars 
represent calculated values of R from experimentally 
measured density of damage. The solid line stands for 
the least square approximation of R^ dependency on J. 
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Resistance Moment Rj'10 (^ 2 ) 

FIGURE 22 

Critical Energy Release rate vs. Resistance Moment 
Solid line represents the average 
value of Gc. 



CHAPTER V 


CONCLUSIONS 


1. We confirmed the history dependence of G c as well as the variance 

of G for the considered loading history, 
c 

2. The experimental observations are in favor of the representation of 

G as a product of the specific energy of damage (which is a material con- 
c 

stant according to CL theory) and the integral cross-section damage 
(history-dependent parameter) 

G - R.Y* 
c 1 

3. Determination of the major energy sink within the hierarchy of 
defects is the critical step for the evaluation of the integral cross-sec- 
tion damage R . According to our estimation, dislocations can not account 

for observed G values. For AISI 304 stainless steel, the slip band extru- 
c 

sions and intrusions are the candidates for being the major energy sink in 
plane stress fatigue. 

4. Our estimation of the value of specific energy appear to be in 
the range of the heat of fusion, in accordance with the discussed 
similarity between fracture and melting (and/or dissociation) processes. 

5. The shape of the crack layer (Figure 16) is similar to the plastic zone in the 
plasticity is helpful in the analysis of stress and strain and, consequent- 
ly, for the calculation of J and M integrals. However, this is only part 

of the problem. The resistance moments R, and R cannot be calculated on 

1 o 

the basis of plasticity. 

6. The crack arrest phenomenon is related to the active zone shape 
change. This transition calls for the definition of constitutive equations 
for (rate of extension) and A (rate of expansion). 
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